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An analysis of the dynamic behavior of a single-pan mechanical balance is presented. In particular, errors 
caused by a swinging pan are analyzed in detail. Results point to a large effect which, though apparently not 
previously appreciated, is nevertheless easily verified experimentally. It is suggested that this effect can be reduced 
to insignificance in a balance whose beam is servo-controlled to an angle perpendicular to the local gravitational 
field. 
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1. Introduction 

The single-pan analytical balance has become standard equipment in laboratories performing precise mass 
determinations. The mechanical balances with which this paper is concerned include those having built-in 
dial weights, a damped beam, and an optical scale which is directly related to beam angle and which has a 
range equal to or- greater than the increment of the smallest dial. Many commercial models of this type of 
balance achieve a precision of better than 0.5 X 10~ 6 of maximum load. Indeed, the most precise kilogram 
comparator — an undamped version of these balances — achieves a relative precision of 1 X 10 -9 [l]. 1 
Equations derived below also treat the case of undamped balances. 

It is a common experience of users of one-pan analytical balances that, occasionally, a badly-centered 
load will cause the pan to swing which, in turn, produces an oscillation of the optical-scale indication. When 
this occurs, one generally discards the reading and repeats the operation with more care. The question we 
posed, however, was: By what physical process does the pan-swing manifest itself as an oscillation of the 
balance beam? The answer might lead to the minimization of the unwanted effect. 

The only previous attempt, to our knowledge, at addressing this question has been in a paper by Bowman 
and Macurdy [2]. These authors reject attempting an analytical solution to the dynamic errors associated with 
a swinging pan because such a solution would be difficult except under unrealistic conditions. Instead, they 
point out that a significant part of the dynamic swing-error is due to the centripetal acceleration of the pan. 
The maximum amplitude of this acceleration is g^ 2 for small ^,, where g is the local acceleration of gravity 
and "j[ is the peak angular amplitude of the pan swing [2]. Even qualitatively, however, this component may 
be excluded as dominant, at least on many analytical balances. The reason is that a swinging pan is seen to 
produce an oscillation in the optical indication at the same frequency as the pan oscillation, whereas the 
centripetal acceleration is modulated at twice the frequency of the pan (i.e., the centripetal acceleration passes 
through two maxima during one complete period of pan swing). 

Our approach to obtaining a solution to the question posed above has been to generate a complete analytical 
solution to the equations of motion of an idealized single-pan balance. Our idealized balance operates with 
knives and flats acting as pivot bearings. The knives do not have perfectly sharp edges but instead have 
finite radii of curvature. This complexity was added to help elucidate the role of imperfect bearings in the 
dynamic behavior of the balance. 

While our balance model, as will be seen below, requires definite restrictions (such as small-angle motions) 
and contains some unrealistic assumptions (such as perfectly parallel knives), it nonetheless includes many 
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of the fundamental physical aspects of the problem and is amenable to analysis. Semiquantitative predictions 
of the dynamic behavior of analytical balances can be made based on solutions of the model. These predictions 
have been verified experimentally as described in section 5 below. 

2. Equations of Motion 
2.1 Method 

The equations of motion for a damped, single-pan balance will be obtained by using the Principle of Virtual 
Work and D'Alembert's Principle [3, 4]. The suspension system, the balance beam, and the counterweight 
are initially represented by point masses located at their respective centers of gravity; they will later be treated 
as distributed masses. 

The first step is to determine the coordinates of these point masses as a function of time. This information 
is used in the Principle of Virtual Work to obtain the generalized torques — and thus the balance beam stiffness 
and sensitivity expressions. The kinetic energy of the system is also obtained from the equations describing 
the locations of the point masses. The kinetic energy, the generalized torque, and a velocity-dependent 
Rayleigh dissipation function [3] are then used in accordance with D'Alembert's Principle to generate the 
equations of motion for the system. 

2.2 Coordinate System 

Figure 1 shows the coordinate system for the balance. Note that the z axis, which is defined to be parallel 
to the local gravitational field, arbitrarily points downward. The y axis is thus parallel to the direction that 
the plane, defined by the two knife edge contact points (P) and (C), would have if this plane were at gravitational 



&— *Y 



(gravitational horizon) 



(B) 




FIGURE 1. The coordinate system and quantities used in the equations of 
motion. The balance beam rotates in the j-z plane, while the suspension 
system-load mass combination, represented by a simple pendulum, swings 
either in the x-z or the y-z plane via crossed knives at pivot point (P). The 
origin is fixed at point (C), which is the center knife contact point for p r = 0. 
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horizon; i.e., if P y = 0. The x axis is directed into the figure so that the coordinate system will be right- 
handed, and the origin of the balance coordinate system is fixed at point (C), which is the center knife contact 
point for P r = 0. Therefore, the balance beam swings in the y-z plane about gravitational horizon. The time 
development of P r (t) describes this motion. 

Points 1, 2, and 3 of figure 1 represent the centers of mass of the suspension system — load mass combination 
{M A ), the balance beam mass (M b ), and the counterweight mass (M B ), respectively. A velocity-dependent 
damping force, F a , is applied at point 4. The knife shown at point (P) is actually assumed to be two knives 
crossed in the same plane (see figure 5 of reference [1] or figure 3 of reference [5]). These crossed knives 
consist of a load knife, L, with a symmetry axis in the x direction and an intermediate knife, /, pointing in 
the y direction. Therefore, the suspension system is free to pivot in any direction about point (P) with an 
instantaneous angle 7. If the knives do not cross in the same plane, then the center of mass length, L cm , can 
be replaced by (i cm )» and (L cm ) y . (The purpose of the intermediate knife in balance design is to compensate 
for non-parallelism of the two main knives [5]. In our model the main knives are assumed parallel, however. 
Thus, in our model the only effect of the intermediate knife is to allow a second degree of freedom for a 
swinging pan.) 

Normally, the suspension system would include a gimbal-like device to minimize the positional effects of 
weights on the scale pan [2]. This would involve a double pendulum on a balance beam — making the problem 
analytically intractable, and also masking the physical effects. (Actually, the gimbal often has enough friction 
to quickly dampen the suspension system motion to that of a single pendulum, so the results of this paper 
still apply to many single-pan balances.) 

We will assume the motion of the suspension system to be that of a simple pendulum. In order to gain the 
maximum physical insight, we consider its motion either to be entirely in the x-z plane or entirely in the y- 
z plane. 

SI units will be used in the derivations, with masses in kilograms, lengths in meters, forces in Newtons, 
torques in Newton-meters, and angles in radians. Note that the quantities y x , 7,, P r , ( b , ( B , and i£ B have 
signs associated with them — and that they are arbitrarily given positive signs in the instantaneous configuration 
at time t shown in figure 1. Small angle approximations: cos"y = 1, cosP r = I, sin", = 7, and sinP r = p r 
will be made at the appropriate places. The 7 approximations are made because precision balances cannot 
function with large suspension system amplitudes, while the @ r approximations are made in order to obtain 
analytical solutions. 

2.3 Coordinates of the Balance 

2.3.1 y m and 0, Mortons 

Assume that the center, load and intermediate knife edges are symmetrically honed cylinders such that 
they have radii of curvature p c , p L and p,, respectively. Then, with the aid of figures 2-4, the coordinates 
of points 1, 2, and 3 of figure 1 are: 

*i = icnSiriT, - p,7, 

y, = - L A cos^ r - P/.sinp T + p c P r 

z, = -L^sinP, - pi(l - cosp r ) - Pc(l - cosP r ) - p^l - cos7.) + L cm cos-i x 

x 2 = 

y 2 = -£ 4 sinP„ + p c P r 

z 2 = £ fc cosP T - pc(l - cosP r ) 

x 3 = 

73 = iflCOsP, ~ ^asinPr + PcP r 

z 3 = Lesinp, + ^aCosP, - pc(l - cosP„), 

where 7, and p r are functions of time. 
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FIGURE 2. Effect on balance beam of center knife rolling on its flat. The 
origin is fixed at point (C), with the positive x axis directed into the figure. 



-j U P sin/3 y 




FIGURE 3. Effect on balance beam of the load knife flat rolling on the load 
knife edge. 
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FIGURE 4. Effect of the intermediate knife rolling on a flat. The positive y 
direction, pointing towards the center knife, is out of the figure. For sim- 
plicity, the edges of the crossed knives are assumed to lie in a plane. 

2.3.2 7 r and P r Motions 

The coordinates of points 1, 2, and 3 of figure 1 are: 

*, = 

V\ = -A4C0S& - PiSin£ r - p t sin7 r + p c P„ + L cm sin7,, 

z \ = -4tsinP y - p/.cos7 r + P/.cosP„ -p^l - cos0„) + L cm cosy y 

x 2 = 

y 2 = -€ A sinp r + p c P„ 

H = VosP,. - pc(l - cosP,) 

H = 

y 3 = I B cosp r - € s sinP r + p c P r 

z 3 = L B sinP r -f- €flCosP r - pc(l - cos P r ). 

2.4 Principle of Virtual Work 

The Principle of Virtual Work [3,4] states that 

W = 2 F,. • df ; = t [MM + (F,)M + (FdM 



where x, = x/7,P r ), y, — y,i"f,fy) and Zj = z,-(7,p T ), and i represents any of the three points shown in figure 
1. 
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Since (F,) x = (FJ r = for this problem, the expression for the virtual work reduces to 



3 



= 2 



™fe) «■ + ™ (t) ^ + ^ (t) * 



(1) 



where Q , Q y and Q^ are the generalized torques for the y x , y y , and (i v motions. It follows from eq (1) and 

figure 1 that 

dz, dz? dz-, 

Qa, = M.g — + M b g — + M B g — . (2) 

Analogous expressions are readily obtained for Q^ x and (^, but we will only be concerned with the conse- 
quences of torques about the center knife, rather than about the load or intermediate knives, because we are 
interested in investigating causes of read-out fluctuations of balances. 

2.5 Generalized Torques 

Using the coordinates listed in section 2.3.1 for y, and (5, motions, and small angle approximations, eq 
(2) becomes 

<? P , = WbbLb ~ MM - [M.APl + Pc) + M^t, + Pc ) + Af rf (€ a + Pc )]p r . (3) 

This expression for Q^ also applies exactly for y r and 3, motions. 

2.6 Balance Beam Stiffness and Sensitivity 
Assume the balance beam is in static equilibrium; therefore, Q^ = 0, and (3) becomes 

[Magi.,, - MM = [M^(9l + Pc) + Mg(( b + p c ) + M„g(l B + p c )]p,, 

which is of the form Tp, = •<£$„, where T^ is the torque about the center knife, ths symbol "*" represents 
equilibrium values, and K is the balance-beam stiffness: 

k = [M A g{p L + Pc) + MM + Pc) + MM + Pc))- (4) 

This result is the single-pan analogue to the stiffness of a two-pan balance with knives of finite radius [6, 

7]- 

The balance sensitivity, S, to changes in torque is 



and has the dimensions of rad/N - m. 2 



s '*t'«- ,5) 



2 Note that if the balance has ideal bearings (p t = pc — 0) and if the position of the counter-poise is adjusted so that (g — 0, the sensitivity is 
completely independent of M 4 an<! *Hj. This result seems to be in contradiction to eq (7) of (7} and to one statement made in that paper. The formulation 
of f 7j can be reconciled with ov.r own if one realizes that the coordinate system of [~\ is itself a function of M E . Hence* for balances whose construction is 
close to ideal (i.e.. t a « {MJM B )t h \, the static properties Isuch as sensitivity) do not depend strongly on load. The dynamic properties (such as period 
of the balance) do depend strongly on load: If a t-kg single-pan mass comparator is converted to a 2-kg comparator by doubling the mass of the counter- 
poise, the sensitivity of the balance will not change appreciably provided the beam and knives do not deform under this modification. The period of the 
undamped balance will increase markedly, however. (A clear and succinct discussion of this point can be found in section 3. A of reference [12]). 
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Note that the balance sensitivity depends upon both € 6 and i B — and that in general neither their individual 
magnitudes nor signs are known. This does not matter, however, since the sensitivity is a directly measurable 
quantity. 

Equation (3) can be rewritten as Q$ r = [MggL B — M^gL^] — K(3 r . Therefore, the effective changes in 
lever arms L A and L B — due to the knife edges rolling on their flats — enters into the quantity KP r ; also, K and 
S will not be constant if the effective values of p c and p L vary with P T . 

2.7 Kinetic Energies 

2.7.1 y. and P, Motions 
The kinetic energy of the system is 

T = 7\ + T 2 + T 3 

T = \m a (x\ + y\ + z?) + \M b {kl + % + %)+ l -M B (x\ + y§ + if). (6) 

Using the coordinates listed in section 2.3.1 in eq (6), employing small angle approximations, and keeping 
only the significant terms, yields the result 

T = \M A [L 2 3 r Y + LL(y,) 2 - 2I cm p,(7,) 2 

+ 2iA„/Y/yA] + \mM,Y (7) 

2.7.2 7, and r Motions 

When the coordinates listed in section 2.3.2 are used in eq (6), the kinetic energy expression becomes 

+ 2L A L cm ( ly + f3,)7 r P r + 2L cm ( Pc - 0J7A] (8) 

2.8 Rayleigh's Dissipation Function 

We will assume that if the balance beam is mechanically or electronically damped, that the damping force 
is proportional to the velocity. Let the damping occur at point 4 of figure 1. Then for small angles, the 
Rayleigh dissipation function [3], F, is 

F = l -a4 = ±aX B & y y , (9) 



and the damping force, F a , is 



dF 

Fa = " *T = ~°^ fl ^' (10) 
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2.9 D'Alembert's Principle 

2.9.1 y, and p r Motions 

D'Alemberts' Pnnclple [3,4] for the equation of motion of the balance beam about the center knife may 
be written: 

d ( BT\ BT BF 

Substituting the expressions for Q^, T and F given by eqs (3), (4), (7), and (9) into (11), and replacing the 
point masses of figure 1 with distributed masses, yields the result 

\l c % + K]0 r + k|J, - [Msgl, - M^ A ] = -[M A L A L cm ](y,) 2 - [M A L A L cm )y x <j x . (12) 

I c is the moment of inertia of the balance about point C of figure 1 when 3 r = 0; i.e., 

Ic = h + h + h (13) 

Ic = [MJft + [(O™ + M&] + [(I B ) cm + M a {L% + €2)] , 

where (/ t ) cm and (/ B ) cm are the respective moments of inertia of the balance beam and the counterweight about 
the x axis through the center of mass, points 2 and 3 of figure 1. The M b (f and M B (L% + €|) terms result 
from the parallel axis theorem. The combination of suspension system plus load acts as a point mass — as far 
as the center knife is concerned — because of the crossed knives at point P. 
Equation (12) has the form 

6p r + 4 r + K0, - d = -eft,) 2 - n,%- (14) 

Appendix A lists all algebraic substitutions used in the derivations, such as those for the coefficients b, c, 
k, d and e of (14). 

The r.h.s. terms of (12) and (14) can be explained by considering the torques about the center knife due 
to the radial and tangential accelerations of the suspension system-load mass combination about point (P ) 
as shown in figure 5. The torque T' c is 

T' c = (^a R cos-y,)(-L A cosP r - p!.sinP r - p c T ) (15 ) 

+ (M^a 7 sin7 : ,)(-L^cosp r - PiSinP, - p c 3 r ) , 

where a R = £ cm ('Y,) 2 and a T = Z. cn ,7,. Note that the torque Y' c is about the pivot point (0) of the center knife 
rather than about the origin (Q, or about the contact point (C); (see fig. 2). For small angles, eq (15) reduces 
to 

T c = ~[M A L A L cm }{y x y - [M A L A L cm ]y x 1 x = -e(-y,) 2 - <rtf„ 

thus accounting for the r.h.s. terms of (12) and (14) by using Newtonian mechanics. 

2.9.2 7, and P, Motions 

Substituting the expressions for Q^, T and F given by eqs (3), (4), (8) and (9) into (11), and using distributed 
masses for the balance, yields the result 
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= -W A L A L cm ](y y )* - [M A L A L cm ]( ly + $ y ) % - [M A L cm (p c - p L )] % , (16) 
which has the form 



b$ r + <rp, + Kp r - d = -e(7 r ) 2 - e( 7 , + P y )7, - J7, 



(17) 



Referring to figure 5, the torques about the pivot point (0) of the knife due to the radial and tangential 
accelerations of the suspension system-load mass combination are 

T c = (^a l ^os7 T )(-£ J< cosp r -p t sinp r - p^inf, - p c r ) 

+ (M A a„smy y )[-L A sin$ y - p^cos^,, + p t co S p„ - Pc(l - cosP,) - p c ] (lg) 

+ (M^ajCos-y^-LrfSinP,, - p L cos7 r + p t cosP T - p c (l - cosP T ) - p c ] 

+ (M^a^in7 T )(-L 4 cosp r - p t sinp r - p t sin7 r - p c p y ) - (M A a T ){-p L ) , 

where a R = L^J^^ 2 and a T = L^fty. The last term in (18) results from a torque about the center of curvature 
of the load knife due to an action-reaction force M A a T applied to the non-slipping contact point (P'); this 
torque is transmitted to point (0) because the balance beam is rigid. 




FIGURE 5. Dynamically induced forces on the swinging pendulum. These forces produce torques about the pivot 
axis of the center knife. This axis, which is in the x direction, passes through the center of curvature of the center 
knife at point (0). 

For small angles, eq (18) becomes 

= -<<7,) 2 -e(7 r +P,)? r -/y r . 

thus accounting for the r.h.s. terms of (17). 
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The — \M A L A L cn \$f) r = — ePjty, term of (16) (which is due to M A times the horizontal component of the 
tangential acceleration times the change in height from gravitational horizon of the load end of the balance 
beam) will later be seen to have a special significance. We will refer to this quantity (divided by I c ) as the 
Mathieu term. (There is also a Mathieu term contribution due to the horizontal component of the radial 
acceleration; but it is negligible for small angles.) The quantity [M A L cm (p c — p L )/I c Yi r will later be referred to 
as the knife-edges term. 

2.10 Physical Pendulum 
2.10.1 Y. and r Motions 

We assume that the motion of the suspension system-load combination is that of a physical pendulum, and 
arbitrarily choose "f x to be a minimum at time t = and to have an amplitude 'f x ; i.e., 

1. = —?,cos(0 7i t , (19) 

where 



< = t ,, , (20) 

'p 



___ [M A gLjj 

I'p = <£). + M A L* m , (21) 



and 



2ir 
w 7, = J~- (22) 

' Tf. 

I'p is the moment of inertia of the physical pendulum about the y axis at point (P) of figure 1, (/s) cm is the 
moment of inertia of the suspension system-load mass combination about the y axis through the center of 
mass point 1 of figure 1, and Py x is the pendulum period. 

The higher order modulation terms have been ignored in (19) because they do not significantly affect the 
balance beam motion. They could be readily obtained by considering the equations of motion about point (P) 
of figure 1. 



Analogous to section 2.10.1: 



and 



2.10.2 y, and P, Motions 



1 r = — ? r c os*V> (23) 



If 



lr = Cs)™ + MaLI,, (25) 



2tr 
w * = jr- (26) 



l„ and (/ S ) OT are the moments of inertia about the x axis through points (P) and 1, respectively of figure 1. 
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3. Solution for y x and ($ y Motions 

The equation of motion can now be expressed in final form by substituting eq (19) into (14), and by using 
the identity cos 2 4i-sin 2 «|i = cos2i|»; it is 

P, + ~ b p r + - b P r - - b = - b -??<c°* 2 <V- (27) 

The general solution of (27) will consist of the solutions to the homogeneous equation (where the r.h.s. of 
(27) is zero), plus a particular solution. 

The homogeneous solution to this equation is easily obtained by standard techniques. The particular solution 
can be found by writing (27) in the form 

d 
P r + u0 y + i?$ r - 7 = wcoskt 
b 

and by assuming that the solution is 

P r = Dcoskt + Esinkt + P r ; 



which yields the results 



and 



uwk 

D = 



[(t^ - k 2 ) 2 + IPk 2 ] 



utf - k 2 ) 



[(t^ - k 2 ) 2 + J?k 2 } 

The general solution of (27) is therefore 

rt 

P r = {« '** [A expVwgt + B exp-V-Wp 2 ,t] 

ct 

or {A + Bt)e~S or 

rt 

e _ 2* [/IcosWp, t + BsinWpt]} 



e , , ( u a - 4w* )cos2&)_ t 

+ e* + i « r — / l < 28 > 



„ ce sin2a) t 

+ 2 - -JX a_ 



[«-4a,^ + 4^<] 



where 



w l = k/6, (29) 

k c 2 



«•>£ = - 



* = I " 47* ' (30) 
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and 

a = d = [MegLs - M&LA 

Pr k [Af^p, + p c ) + M„g(t b + p c ) + M^(€ fl + p c )] 

The choice of terms in eq (28) arises because the balance may be overdamped, critically damped, or 
underdamped depending on whether <?/4b is greater than, equal to, or less than k. For all damping conditions 
the motion of the balance beam about the equilibrium angle |3 r is modulated at an angular frequency which 
is twice that of the natural angular frequency ta yK of the physical pendulum. This modulation effect is due to 
the radial and tangential accelerations of the suspension system-load mass combination. 



4. Solution for y y and |J r Motions 

4.1 Equation of Motion 

The equation of motion is obtained by substituting eq (23) into (17), and by using the identity cos 2 «ji — 
sin 2 !)* = cos2»|*: 



Wz + i 
d? b 



d$ y (k e \ d e f 

~27 + \b + l % * coso vJ p ' = l + l ^ W > cos2w ^ ~ 6 %0i *- cos<v " {32) 

The I? V *, cosw 7y t jp T term of eq (32) — which is the Mathieu term referred to in section 2.9.2 — complicates 

the problem; yet it cannot, in general, be neglected when considering the dynamical effects of precision 
balances. Therefore, (32) will be converted to a more tractable form. 

First, we change from the time variable, t, to the dimensionless variable z by making the substitution 

<V = 2z - (33) 

Equation (32) thus becomes 



*• +2 






4{% 



cosu, t , (34) 



or 



£0, -rfp 

dz 2 + 2k ~d7 + (° ~ 2 9 cos2 *)Pr = h + ^cos4z - mcos2z = ~f{z). (35) 

We then replace P r with the quantity [10] 

r = e- b u{z) ; (36) 

thus obtaining the equation of motion in final form: 

d 2 u t 

j-J + (« - 20cos2z)u = e h f(z), (37) 
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where 



= a - k 2 



_L \l _ ill 



(38) 



The quantity a is negative if the balance beam is overdamped, zero if critically damped, and positive if 
underdamped; whereas q is always a negative number. The quantities z, f(z), u(z), k, a, a, q, h, (, and m 
are all dimensionless. The general solution of (37) consists of the solutions to the homogeneous equation, 
plus a particular solution. 



The homogeneous equation 



4.2 Homogeneous Solutions 



u" + (a — 2qcos2z)u = 



(39) 



has the canonical form of Mathieu's equation [9,10]. The reader is referred to reference [10] for an excellent 
account of how to obtain solutions to this equation. 

Figure 6 shows an isostability diagram [10] of v for possible solutions to eq (39). We are only interested 
in underdamped or critically damped solutions because laboratory analytical balances are seldom designed 
to operate in the overdamped mode; therefore 1 >> a > 0. Also, the quantity q = —2M A L A L cm "i y ll c is 
always a very small negative real number for precision balances; so — 1 << q < 0. The region of interest 
for these values of a and q is indicated by an arrow in figure 6. Clearly, the solutions of (39) cannot be 
"simple" Mathieu functions of the sine or cosine types s(e^), c(e ) or c(ej) [9,10], but instead are more 
complicated functions. 



c(ej 



c(e 2 ) 




FIGURE 6. An isostabilily diagram of V for possible solutions of Mathieu's 
equation given by (391. The arrow points to the region of interest for precision 
balances; i.e. 1 » a & and — 1 << q < 0. Solutions in the shaded 
region are unstable. This does not mean that an overdamped balance (a < 0) 
is mechanically unstable because u(r) is multiplied by the factor e" 4 " in the 
solution for P,. 
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Since q is small, we guess that the solution is of the Floquet form [9] 






where v is not an integer; (i.e., u(z) is not a Mathieu function). For small values of a, v is defined by [9] 

/ aw 2 a 2 ir 4 \ fit 2 [ / tt 2 \ 1 

cosvlr = ^ _ _ + _ + . . . j _ £_ ^ + a ^ _ _j + . . . J 

V \96 256 / 

We then use the identity e' 9 = cos0 + i sin8 in (40) to obtain the two independent solutions of eq (39); 
they are 

u x = cosvz — qcosVzcos2z — v<7sinvzsin2z (41) 

and 

u 2 = sinvz — gsinvzcos2z + Vgcosvzsin2z, (42) 

where terms of order f have been neglected, and 

" 9 = 27I^v < 43) 

v « \/^ (for a > 0), (44) 



v~^ (fora = 0). 

The homogeneous solution of (39) is therefore 

A Ul (z) + Bu 2 (z). (45) 

4.3 General Solution 
4.3.1 Particular Solution 
We must next find a particular solution, u p {z), of 

u + (a - 2 9 cos2z)u = e**/(z), (46) 

which has the form of Hill's equation [10]. The solution is obtained from the expression 

",(*) = ~ ^ |^,(z) j o u^z^ftfdz - ujz) j' o Ui(*)e^)A 1 (47) 

where u, and u, are given by eqs (41) and (42), /(z) = h + €cos4z - mcos2z, and 
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? = Ul (z)i4(z) - u,(zK(*) = v. (48) 

The evaluation of eq (47) is straightforward, but very lengthy; the result is included in the next subsection. 

4.3.2 Complete Solution 
The general solution to eq (35), using (36), (45), (47), and (48), has the form 

P r = e-' h u(z) = e-Hlu.W + Bu^z)] - '-^ «,(*) §u&Wfc)dz - u 2 (z) £u 1 (z)e%) £ £ Z J, (49) 

which, when using (41) and (42) is 

P y = i4e -fa [cosvz — qcosvzcos2z — vgsinvzsin2z] 

+ fie~**[sinw? — 9sinvzcos2z + v^cosvzsin2z] 

- w^i{{ h + 2 m ) I_v + e ~* (vcosV2 + ~ Mnvz)] 

— (q/i)[ — vcos2z + e -fa (vcosvzcos2z + feinvzcos2z)] 

— (vg/i)[fcin2z + e _ * s (vsinvzsin2z + Acbsvzsin2z)] > 

= L_ J ( _ m _ a h - 1 ( _ V g h + Yl ( )[-fein2z 

- (v-2)cos2z + e" fe ((v-2)cosKz + Asinvz)] 

+ (9"i)[-fain2zcos2z - (v-2)cos 2 2z + e _fa ((v-2)cosvzcos2r + feinvzcos2z)] 

+ (v9m)[-(v-2)sin 2 2z + Jfcsin2zcos2z + e _fa ((v-2)sinvzsin2z - *cosvzsin2z)] \ 

= 1 { ( -m - qh - 1 1 + Vqh - V J € ) [*sin2z 

2v[F + (v+2) 2 ] H q 2 H 2 / 

- (v+2)cos2z + e- b ((v + 2)cosvz + fcsinvz)] 

+ (gm)[Lin2zcos2z - (v+2)eos 2 2z + <r fa ((v + 2)cosvzcos2z + fcinvzcos2z)] 

+ (vgm)[Asin2zcos2z + (v + 2)sin 2 2z + e -**((v + 2)sinvzsin2z - Acosvzsin2z)]J 

= 1 \(i + £ m + v?mj[-fain4z - (v-4)cos4z 

2v(fc 2 + (v-4) 2 ] \\ 2 2 / l 

+ e~ fa ((v-4)cosvz + fainvr)] 

- (?€)[- Jfccos2zsin4r - (v-4)cos2zcos4z + e- 6 ((v-4)cosv*cos2z + fainvzcos2z)] 

- (vg£)[fein2zcos4z - (v-4)sin2rain4z + *- fa ((v-4)sinvzsin2z - Icosv»in2z)]| 

= I l(( + i m - v * m ) [fein4z - (v+4)cos4z 

2^ + (v + 4) 2 ]\V 2 2 / l 

+ e~ fa ((v + 4)cosvz + fcsinvz)] 
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- (g€)[fccos2zsin4r - (v + 4)cob2zcob4« + <.-*((v + 4)cosvzcos2z + *sim*cos2z)] 

- (vg€)[fain2rcos4z + (v + 4)sin2zsin4z + e~ fa ((v + 4)sinvzsin2z - Acosv*sin2z)] > 

7= ~, ^ I I - ^ - v I < )[-fain6s - (v-6)cos6z 

2v[fc* + (v-6) 2 ] \\ 2 2 ) V 

+ e _fa ((v — 6)cosvz + feinvz)] > 

7=-, — -, Wi\ ( - - € + v 2 ( )[ksm6z - (v + 6)cos6z 

2v[fe 2 + (v + 6) 2 ] 1\ 2 1 J 1 

+ e -£l ((v + 6)cosvz + fcsinvz)]>. (50) 

This is the exact solution (in reduced notation) for small angle 7 and @ r motions. It is obviously much more 
complicated than the solution for small angle "y^ and P r motions given by eq (28). 

It requires too many printed pages to provide an expanded form of eq (50) written in the notation of (28). 
(There are 74 terms in that notation.) However, the reader can readily transform any. term of (50) into that 
notation with the aid of Appendices A and B. Equation (50) is provided so that readers seeking solutions to 
problems similar to our own may be spared the considerable tedium of its derivation and so that the terms 
neglected in the sequel may be recovered by the interested reader. 

5. Special Cases 
5.1 A Freely-Swinging Balance 

In a completely undamped balance, k = c = 0, and Wp T = <Op v ; so eq (28) becomes 

P T = 4cosu) Pr « + BsinWpt + p r + - -y 2 2 J" cos2co 7x t. (51) 

Using Appendices A and B, and keeping only the largest terms, (50) becomes 



P r = /icoscop.t + BsinWpt + 



P -^"f- 



2b 2 ' y w|,(^ r ~ 4<4 



>] 



6 ((o 2 r - 4w| v ) 7 " b (w|, - W 2 ,) 

+ " "?r T~, — * — TT cos2u>, t + — H — - — *-r — cos3to, t, (52) 

6 « - 40 * 26* 7 ' « - 4a> 2 ,)« - 9«y 

where 6 = I c , e = M A L A L cm and f = M A L rm {p c — p L ). Only cosine modulations of the /tcosa)p r f + 
Bsino)p f motion appear in (51) and (52) because of the choice of the phases of the y, and *y, motions given 
by eqs (19) and (23), respectively. 

The cos3(0 7r t term of (52) is much smaller than the cos2o) yj t term, and is therefore negligible. There is a 

systematic error term in the equilibrium angle: — e/"^^o>*/[26 2 <»)| r (u)^ — 4u>|,)]. This term is probably small 
for most precision balances and may in fact be an artifact of our initial neglect of terms in "ft. 

The size of the costo> 7r f modulation depends upon the knife edge radii p c and p t , and also upon the value 
of P r The cosw 1r i modulations will be at a minimum for a given "^ y , when P, is chosen (usually via an 
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appropriate adjustment of M A ) such that the Mathieu term contribution is nearly equal and opposite to the 
knife-edges contribution. These two terms then nearly cancel, and one should observe pure cos2a>, r f mod- 
ulations. In most realistic cases, the knife-edges term can be safely ignored. (The radius of curvature of a 
good knife is usually less than 1 u.m [12].) Then the Mathieu term vanishes when the equilibrium position 
of the beam is parallel to the gravitational horizon (P v = 0). Note that a test for the absence of a systematic 
error contribution to (J r does not a fortiori assure a negligible knife-edges contribution to the cosw^t modulations 
because the error term goes as ^ 2 , whereas the modulation term depends upon "V, — and is therefore much 
larger. 

5.2 A Damped Balance 

5.2.1 Theoretical Predictions 

We shall now investigate the motion of a damped balance at a time t which is long enough so tha 

involving the factor e~** = e~ at2h become negligible; (i.e., a time such that, in the absence o 
the balance beam would be at rest at the equilibrium angle 0,). Equation (31) thus becomes 



P r = P r + ~ h t 



Ice 



<« - 4Q 



(io| r - 4(0 2 ,) 2 + 4 



H 



cos2d) Yi t 



+ — T 



w" 



(< - 4< ) 2 + 



4c 2 1 



sin2to > t, 



(53) 



and (50) becomes 



K = [P. - 



« -«. 



r. 



2b 2 * <4,« - 4<o£). 

2 



io; 



+ b P ' % « - 4<) 



- 1 r •= ^ =: cosw 7 ,t 



_0i_ 



&2 Pv1r (< - 40$) 



b 2 



12 % 



O 2 + 



N 



+ _*2_ 
b 1 " r 



<K, - 4<) 



« - 4<) 2 + 4 



U] 



cos2(i) y t 



Ice 
I 1 



+ -? 2 



sin2w v t 



(w 2 r - 4w 2 ,) 2 + 4 



$<] " 
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«j2 G>! ((1)S - 9o£ ) 

+ £ 15 F*^ ** 7"! cos3a V 

«-4o,^«-9<) 2 + 9-<J 



See 2 ... a> 5 
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+ — f y p * =■ S in3o) 7r t, (54) 



« - 4<)|\< - 9<) 2 + 9 ^ <j 



where b = I c , c = a^l, e = M A L A L cm ,f = Af^Jpc-pJ, and Wp r = at critical damping. 

Sine modulation terms now appear in (53) and (54) because the damping effects introduce phase shifts into 

the balance beam responses to the pure cosine y, and "y T motions. The sin — w^t terms are smaller than their 

n ■ . -2c 

equivalent cos — ti)~t terms by the ratio — - , where n — 2, 4, or 6. 

2 """V 

The 3o Tr t modulation terms are negligibly small compared with the 2<a y t terms — as was the case for a 

freely-swinging, undamped balance. The presence of a systematic error term in the equilibrium angle can, 

in principle, be tested by observing the apparent P values for pure y, and y r motions. In practice, however, 

the over-simplifications of our model may cause problems in interpreting such a test. The modulations of 

oi . 

frequency — " can again be minimized by choosing a nominal value of P T such that terms involving the 
2ir 

quantity M A L A L cm $ y 'ii r are nearly equal in magnitude and opposite in sign to the knife-edges terms which 
include the factor M A L K Jp c - pjy,. 

5.2.2 Experimental Verification 

We have tested some of the predictions of section 5.2.1 by using a Mettler M5 balance, 3 which has a 
velocity-dependent damping produced by an air-dashpot. This balance has a maximum load of 20g, an on- 
scale range of 20 mg, and a readability of 2 pg. An analysis of some important features of this balance can 
be found in ref. [7], The weighing pan is attached to the lower part of the suspension system by a hook and 
ring, so that the initial motion of the suspension system is that of a double pendulum. This motion quickly 
dampens to that of a single pendulum — which is a necessary condition for our model. The presence of the 
hook and ring, however, makes it difficult to generate pure 7, or pure y y motion. We therefore temporarily 
locked the hook and ring together to provide a rigid suspension system for these tests. The tests described 
below were performed by gently starting the pan swinging at an amplitude of 4 ± 1 mrad about equilibrium. 

What can be inferred from eq (54) about the dynamic behavior of this balance? First, eq (54) can be 
simplified because <o^ = and w^ « a> 2 : 



P r = const. 



+ 






-;![ i + fe)H 2 "*' + Area "^) (C) 

+ higher order terms. (•*>) 

Brand name* are used only for purposes of identification. Such use implies neither endorsement by the National Bureau of Standards nor assurance 
thai the equipment is the best arailablc. 
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Term A, the Mathieu term, arises from the tangential acceleration of the swinging pan. Term B is a consequence 
of having knife edges of finite radius. The last term considered, C, has its origin in the centripetal acceleration 
of the pan. The latter is the only term of importance in eq (52). 

We have made the following estimates for the balance being considered: 

gem 2 

|/| *£ 2.5 X 10- 2 gcm 2 
b ~ 3.5 X 10 3 g cm 2 
c~ 5.5 X 10 3 gcm 2 s- 1 . 

Also, u>^ was measured to be 7.8 rad/s. 

On the basis of these numbers, one may observe that term B always has much smaller amplitude than term 
C and is, therefore, always masked. We are left only with terms A and C. Equation (55) predicts that term 
A will dominate when |p y | » y y but that term C will dominate when the inequality is reversed. Thus, under 
the condition that |P y | >> y y we would expect an oscillation in the beam at the same frequency as the pan 
oscillation whereas the oscillation frequency of the beam should double when |3 y = 0. The ratio of the amplitudes 
in the two regimes is predicted to be 

4.4 $% . (56) 

Note that this result is independent of the estimate of e given above and only weakly dependent on c and 6. 
These predictions were tested in the following way. 

With the balance arrested, a load of nearly 20g was placed on the pan. When the balance was released, 
the optical scale reading was found to be near zero (i.e., 8, was near its maximum value of +22 mrad [7]). 
Pure y x motion produced a small modulation of the balance beam read-out angle, B v , at a frequency twice 
that of the suspension system — as predicted by eq (531 — at an amplitude of 2 ± 1 p.g. Pure y motion 
yielded a modulation having an amplitude of 32 ± 2 u.g — with a frequency which was the same as that of 
the suspension system. The phase shift due to damping effects was too small to observe, so that this modulation 
was a nearly pure cosa> 7 f function. Our measurements could have detected a phase shift as small as Tf/6 but 
eq (55) predicts a phase shift of about tt/16. 

We then changed the sign of 0, by adding 20 mg of tare weights, thereby moving the equilibrium of the 
beam near to the maximum on-scale reading (i.e., (3 V near its minimum value of —22 mrad [7]). The result 
of pure "y v motion at this position produced a 180° phase shift in the cosu) 7 ,t modulation of B v as predicted 
by eq (55). 

The magnitude of M A was then reduced by removing 10 mg so that the equilibrium position of the beam 
was at the middle of the optical scale. At this position, where the beam is nearly horizontal, the Mathieu- 
term contribution of the cosw 7 ,f modulation nearly vanishes. The 8, modulation now became that of the small 
cos2(o^t contribution (in agreement with the case of y, motion.) 

These tests clearly verify several important predictions of eqs (53) and (55). The observed ratio of amplitudes 
in the two frequency regimes is 16, which agrees satisfactorily with the value 24 predicted by (56). 

6. Conclusions 

In summation, it will be useful to recall some of the important simplifications in our mathematical model. 
All knives have been assumed to be symmetrically honed and to roll on perfect flats. The two main knives 
of the balance are assumed to contact their respective flats in straight lines which are parallel to each other 
and to the x-y plane. Many authors have dealt with static errors which result from a violation of those conditions 
(e.g., [5, 11, 12]). An analysis of the corresponding dynamic errors is clearly beyond the scope of this paper. 
We have also chosen to look at pan oscillations in the context of the small-angle approximation of an idealized 
balance. While we believe this is appropriate, it should be mentioned that small variations in balance sensitivity 
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as a function of (3* can be observed in many analytical balances having a large on-scale range [7]. Such 
behavior cannot be predicted from our equations, which ignore angular terms higher than first order. 

The above considerations serve as a caveat to the quantitative application of our results to real balances. 
Nevertheless, the solutions found have proved valuable. The question posed in the introduction has been 
answered: Two processes are important in analyzing the effect of a swinging pan on balance read-out. One, 
which has been mentioned in the past [2], is due to centripetal acceleration of the pan. The second — and 
often predominant — process is a more subtle effect by which the tangential acceleration of the swinging pan 
leads to a torque on the beam if the latter is not at the gravitational horizon and if the former has a component 
of motion in the y-z plane. Even a crude approximation of the tangential acceleration (e.g., a^max.) ~ 
g*f) suggests that its influence will exceed that of the centripetal acceleration if (3 > 7. A convincing argument, 
however, can only be advanced through the rigorous mathematical solution to the equations of motion as 
presented in section 4. 

The disappearance of this second, or "Mathieu" process when the rest-angle of the beam is near gravitational 
horizon may be used to minimize the effect. Alternatively, it may be possible to use the effect to adjust the 
beam angle in a class of experiments [13, 14, 15] where it is advantageous that the beam be parallel to 
gravitational horizon. Here the radii of the knives must be determined (under load) in a separate experiment 
[12]. 
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Appendix A. 
Algebraic substitutions and changes of variables 

This appendix lists all of the algebraic substitutions and changes of variables used in the derivations. 
Therefore, any equation given in the main text can be readily expressed in terms of the basic quantities y x , 
7 r , $,, L cm , L A , €/,, L B , !£ B , ( B , M A , M h , M B , g, a, t, p c , p L , and p, — as defined in figures 1-4. 
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* = Wa g(pL + Pc) + M h g(e b + Pc ) + M B gi£ B + p c )] 
1 



5 = - 

K 



I c = MJ 2 + [(/ 4 ) cm + MJ*JL + [(/ fi ) cm + M^Ll + m 

Ip = &)„ + M A l* 

I'p = (lsl m + M A Ll 

,2 _ i M AgL cm ] 

»y. -—7, 

*p 

W A gL cn \ 



I 2 

cm 

I 2 



0> 2 = 

7r I 



w y. = 



<p 

Py, 

2ir 



ft),, = 

Py 



"1 X — ~ % COS <*yj 

y r = — % cos o> 1v / 

6 = / c 

c = a^I 

<* = [^figifi - M A gL A ] 

e = M^LAm 

/ = MaLcJPc ~ Pi) 
A _ d _ [M B g L B ~ MAgL A ] 

' * [M A gip L + Pc) + Mtfitt + Pc) + M B g{€ B + p c )] 



b 46 2 PL 



Pv 




[(r 2 - F) 2 + U 2 * 2 ] 
£ »fr» - **) 
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_ C 

* = r — = 

0Ci) 7> , 


oSPb 







e . MJLJ., 



= -2^ T = -2^-^r 



g =-z- 7r --. 



c 



h ~ <6 / c < 

O If 



lr 



/(z) = A + £cos4z — mcos2z 
Pr = e"* u(z) 



4 

a = S- *?« — 



6 46 2 J 



*■ 2 



<oj», = (if a = 0) 

- _ g _ _ e (Q^-y v 

9 ~ 2(1 - v 2 ) 6 (< - 4(o'|,) 



v * >A » 2 ^ (if a > 0) 



v~^|=V2g= -V2^% («fa = 0) 

« t = cosvz — g cosvz cos2z — vq sinvz sin2z 
Uj = sinvz — q sinvz cos2z + vq cosvz sin2z 
p = u t i4 — "2 U 1 = v 
vz = w p P ' 
2z = <«V 
I - ct 

g <*>-,,<%?, 



Appendix B 
An aid to identifying terms in eq (50) 

This: appendix provides some relationships that can aid the reader in transforming the terms in eq. (50) 
for y r and (3 7 motions into the notation used in (28) for -y x and {$ r motions. The combination 

(BD 



[¥ + (v + n) 2 ] [fc 2 + (v - n) 2 ] 
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often occurs in (50), where n = 0, 2, 4 or 6. By using the common denominator 

[(5 . „ T _ ^ - ^ [(* - £ <)* + i £ <]. 



(Bl) can be written as 



nvktoi 
* (if X = k and y = -A); (B2) 



as 



{M<)* + £.<] 



jfc(4a)i + nW )c> 2 
- L — §I 2lZ - Jl (if X = k and Y = A); (B3) 



•[(•*- t<)* + tN 



-*— & *£_Jh (if JiT = ( V + n) and y = (v - n)); (B4) 



'[(«*- T<)" + 7?<] 



or as 



n (4o)^ + n 2 < - 8o>'|> 2 



[(+-H ,+ i?<] 



^ (if X = (v + ii) and Y = -(v - n)). (B5) 



With the help of eqs (B1)-(B5) and Appendix A, the terms of (50) can be expressed in the notation of 
(28). For example: 

mk sin2z mk sin2« cf <a% 

= - rz %-z * «n o> >( , 



2v [k 2 + (v + 2) 2 ] 2v[k 2 + (v - 2f) 6* 



[«~<) 2 + ^<] 



and 



vqhk sirOz vqhk sin2z ce P r o)?, w„ (wf, + wi ) 

2v[A* + (v + 2)*] 2v[A* + (v - 2f] h^> |" ,„.«» t l * 
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